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Abstract The editorial handling of papers in scientific journals as a human activity 
process is considered. Using recently proposed approaches of human dynamics theory 
we examine the probability distributions of random variables reflecting the temporal 
characteristics of studied processes. The first part of this paper contains our results 
of analysis of the real data about papers published in scientific journals. The second 
part is devoted to modeling of time-series connected with editorial work. The purpose 
of our work is to present new object that can be studied in terms of human dynamics 
theory and to corroborate the scientometrical application of the results obtained. 
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Introduction 

Nowadays new possibilities for quantitative studying of human activity processes ap- 
pear. Modern computer technologies allow to collect and to process large volumes of 
statistical data. It is possible to fix the time of each performed operation from large flow 
of human (customer/user/client) actions and to obtain the data set representing whole 
picture. Therefore, new methods for data analysis and knowledge discovering could 
be applied. In particular, this allows quantitative analysis of psychological, social and 
other aspects of human behavior. 

The sequences of human actions (telephone calls, information queries or stock ex- 
changes) are not new subject to study. But contrary to generally accepted opinion 
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about human actions randomly distributed in time and thus weU approximated by 
Poisson processes (interevent times are exponentially distributed), the comparatively 
recent results of human dynamics analysis show the presence of power laws which na- 
ture and origin are unexplained up till now [T][21[3]- This discovery attracts an interest 
and provokes analysis of new processes involving human actions in different fields of 
our life. The analysis of time statistics of human activity patterns can be useful for 
different optimization and control tasks in spheres of mass service, communication, 
information technologies, resource distribution, etc. Besides, this approach can give a 
new possibility to understand human behaviour and to get its additional quantitative 
measure. 

The continuous human behaviour can be considered as the set of consecutive actions 
in time. In this case it is convenient to examine two characteristic random variables: the 
time interval between two consecutive actions (called the interevent time tint) s^nd the 
time a task is waiting for an execution (the so-called waiting time tw)- The power-law 
nature of corresponding distributions 

P{tint)-tr^, P{t^)^t-P, a,/3>0 (1) 

arises from the analysis of real processes that involve human actions such as browsing 
the Internet, data downloading, electronic and mail communication, initiating financial 
transactions etc. [T][2l[4]. On the other hand, the processes with execution of tasks in 
queue are the typical objects of study in the queueing theory which can be seen as 
the subfield of the so-called theory of mass service or, generally speaking, of applied 
probability theory [5l[6]. So, it is natural to use some general terminology to describe 
such kinds of processes. An empirical analysis results in a variety of values of the 
exponents governing ([T]) [3115], however an observed tendency when similar values of 
the exponents describe different processes gives rise to an analogy with the "universality 
classes" found in physics of critical phenomena [5]. Different models have been proposed 
to gain insight about an origin of power laws IT]) [3]. Some of the proposed models are 
based on the assumption about the key role of priority-based (decision-based) queuing 
process [5] . The situation when an individual has to perform a list of tasks and chooses 
a task from this list using some internal priority is very natural. It is important, that 
power-law-like distributions ITJ appear only in the so-called critical and supercritical 
regimes of service, when its traffic intensity p: 

p = A/m (2) 

is greater than 1. In (12]) A is the tasks arrival rate and p, is the execution rate, re- 
spectively IS]. In other words, the probability distribution functions of human activity 
processes are close to power law only when the queue of tasks is not exhausted 15]. 
This class of human activity can be described by the so-called "task-driven" models 
[3]. On the other hand, the "interest-driven" models explain the existence of power laws 
in human dynamics in a different way. There, the power laws may be observed also 
for the processes without any possibility to determine the queue of tasks, for example 
visiting of web-sites by different users 3j. 

Thus, we can conclude that human dynamics theory is the issue of current impor- 
tance. Often besides academic interest there is an obvious great practical value of the 
results obtained. 

The main purpose of this paper is to gain insight on the editorial process in sci- 
entific journals using tools and concepts of the human dynamics theory. In the first 
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part (section [T]) we show our empirical results obtained during the data analysis about 
distributions of waiting times of papers in several scientific journals. The simple simu- 
lation model of editorial processing of scientific manuscripts is presented in the second 
part of this paper (section [5}. 



1 Waiting times statistics analysis of scientific journals 

Let us consider the editorial process in scientific journals as an example of human 
activity processes [?]■ Keeping in mind the classification mentioned above, we will 
consider the "task-driven" model. In this kind of mass service system the input flow 
consists of submitted papers forming the queue. A standard procedure after paper 
submission can include the following steps: (i) peer-review process, (ii) revisions if 
necessary, (iii) acceptance by an Editorial Board, (iv) other intermediate processes. On 
each of the above stages the paper may be rejected. However, typically the information 
concerning the rejected papers is not publicly available. Therefore, we consider the 
random variable tw defined as a number of days between the dates of the paper final 
acceptance ta and the paper receiving tr: 



All stages of the editorial processing of submitted manuscripts are considered together 
as the one process (Fig. [T|. Though more than one actor takes part in it, we consider 
that every part of this work is controlled by an Editorial Board. Therefore, we can 
treat this process as one of the main characteristics of the Editorial Board activity. 



received papers 



accepted papers 




date of receiving, f, date of acceptance, 4 

Fig. 1 Schematic picture of editorial processing of manuscripts in scholarly journals. 



During the real data analysis we met the problem of acceptance time ta, determining 
whereas the dates of paper receiving tr are usually fixed very accurately. Different dates 
refiecting stages of papers processing are available for different scientific journals: date 
of revision, date of final acceptance, date of availability on-line, etc. It is necessary to 
specify the meaning of ta for every particular journal. Besides, it is also interesting to 
consider various sets of tw choosing different meanings of ta for the same journal and to 
compare the obtained results. In our analysis, the ta was defined as the date of revised 
version if the final acceptance date was omitted and as the final acceptance date (if 
present). So, tw were calculated as the time differences between the date of submission 
and the most distant of two possible dates: revision date or final acceptance date. 

Our goal was to determine the functional form of probability distributions P(tw) 
based on the statistical data analysis performed for a few scientific journals. Another 
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task was to find if possible a typical form of P (tw) for normally working Editorial 
Board. 

Several journals with different Editorial Boards were chosen in our study: three of 
them belong to the international Elsevier Publishing House ("Physica A: Statistical 
Mechanics and its Applications", "Physica B: Condensed Matter" and "Information 
Systems") [7]- One more journal is rather new "Condensed Matter Physics" (published 
by the Institute for Condensed Matter Physics: www.icmp.lviv.ua). The publicly avail- 
able data from the official web-sites were used for analysis (Fig. Also we tried to 
get the corresponding statistics of "Scientometrics" , although the data set is too small 
because the required dates for papers are available online only for several last years. 
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Abstract 





Fig. 2 An example of publicly available data (the official web-site of the journal "Physica A") 
which were used for our research. 

Some formal parameters of databases used in this part of our work are shown in 
table [1] As we can see from the table, the typical number of records for each journal 
is of order of 10^ which allows to make some quantitative conclusions. These general 
conclusions hold even for a more poor statistics (c.f. results for "Condensed Matter 
Physics" in Fig.|Bl). Zeros in the table mean that for some papers their submission date 
is coincide with the final acceptance date. But what is even more interesting, in every 
journal one can find papers with very long waiting periods (for example, 2260 days it 
is more that 6 years!). Oh course, we can only speculate about the possible reasons for 
that. For example, it could be long discussion about the manuscript or banal misprint. 

At the first stage the probability histograms of waiting times P{tw) for selected 
journals were constructed fjl . All experimental values were distributed among discrete 
intervals (bins) of length equals to 5 days. 

Our purpose was to verify the functional form of P(tw) and to refer it to the power- 
law-like class (non-Poisson processes) or, for example, to the exponential-like class 
(Poisson processes). The exponential distributions of random variables tint and tw are 
evidences of the random selection of tasks to execute [3ill0| . In Fig. [3] we see the P(tw) 
distribution for journal "Physica A". It has a form of a unimodal non-symmetrical 
distribution with a smooth decay. Since the distribution is skew, the mean value f™"*^" 
does not coincide with the typical value fw^, at which maximum in P(fw) occurs. Both 
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Table 1 Characteristics of the data sets analysed: general number of records, maximal (t™'"'), 
minimal (t™'"), typical (t^^), mean (t™"'"^), and median (t™'^'^) waiting times for the journals 
under consideration. 





"Physica A" 


"Physica B" 


"Information Systems" 




(1975-2010) 


(1988-2010) 


(1975-2010) 


Number of records 


4576 


4944 


814 


C"", days 


1629* 


1087* 


2260* 


C". days 











iS'P, days 


60 


80 


245 


C'''"', days 


124 


122.2 


331.7 


C^'', days 


95 


90 


275 



* LA. McLure, A.-M. Williamson, Physica A, 1996, 234, Iss. 1-2, 206-224; ibid. 225-238. 

* V.G. Bar'yakhtar, V.A. Popov, Physica B, 1999, 269, Iss. 2, 123-138. 

t M. Binbasioglu, D. Karagiannis, Information Systems, 2000, 25, Iss. 6—7, 453-463. 



the mean, typical and median values of tw are given in the table \T\ The log-log and 
log-linear plots of Fig. [3] demonstrate good possibilities for linear approximations of 
P{tw) in both scales and this situation is analogous for other journals analysed [7]. 




Fig. 3 (a) Log-log plot and (b) linear-log plot of the P(tw) distribution for journal "Phys- 
ica A" . 



Further, we have verified two main hypotheses about the form of probability dis- 
tributions which are used to describe human activity processes: (i) log-normal distri- 
bution [5]: 

A [-(If)]' 

P{t^) = Po + ^= e , tc,Lu>0, (3) 

where In (tc) and cj are the mean and standard deviations of the In (tw), Po, ^ are fit- 
ting constants; and (ii) power-law distribution with exponential cutoff [5] for exponent 
values a — {l;3/2}: 

P(tw) = At-"e"^, to>0, (4) 

where to is characteristic of waiting time which depends on traffic intensity, A is a 
constant. Using fitting procedure we found optimal parameters for both distributions 
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(O and (|4]). The results of fits obtained for journal "Physica A" are shown in Fig.|4]by 
smooth curves. 




Fig. 4 The P(tw) distribution for journal "Physica A" with different approximation curves: 
(a) log-normal ^ (solid line, red online), (b) power-law with exponential cutoff l(4j with a = 1 
(light solid line, green online) and a = 3/2 (dashed line). 




200 400 600 800 1000 400 800 1200 1600 2000 2400 



Fig. 5 The P(tw) distribution for (a) "Physica B" and (b) "Information Systems" journals. 
The dotted dark (blue online) lines are the approximations by log- normal JSj, the light solid 
(green online) and dark dashed (red online) lines - by power-law with exponential cutoff ||4j 
with a = I and a = 3/2, respectively. 



To compare the accuracy of approximations by different functions (log-normal ([3} 
and power-law with exponential cutoff (Q) we exploited the criterium based on the 
value of the adjusted coefficient of determination [§_. . This coefficient is used to 
verify the closeness of experimental data to the non-linear theoretical curve and is 
a modification of statistical coefficient of determination which is the square of 
the sample correlation coefficient between the outcomes and their predicted values. 
A value of R^ close to 1 indicates that the fit is a good one. For example, for the 
journal "Physica A" we have found that both log-normal ^ {R? ~ 0, 97) and power- 
law function with exponential cutoff @ and exponent a = 1 {R? « 0, 95) can be the 
probable functions of distributions P{tw)- The value of R^ for power-law approximation 
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function with exponential cutoff (|4]) and exponent a = 3/2 is slightly smaller (~ 0, 92) 
for this journal. 

In fact, both log-normal and power-law functions predict the same leading behavior 
t~^, differing only in the functional form of the exponential correction [TD]. The advan- 
tage of hypothesis about log-normal functional form of P(tw) consist in the possibility 
to describe all the data span, not only the tail. The results of analogous approximations 

0,050 -j- 

0,045 - 

0,040 - 

5 0,035- 
Q. 

0,030 - 
0,025 - 
0,020 - 
0,015- 
0,010- 
0,005 - 
0,000 - 
-0,005 - - 



Fig. 6 The P{tw) distribution for "Condensed Matter Physics" journal. The dotted dark lines 
(blue online) are the approximations by log- normal (|3jl, the light solid (green online) and dark 
dashed (red online) lines - by power-law with exponential cutoff (|4]l with a = I and a = 3/2, 
respectively. 

obtained for "Physica B" and "Information Systems" journals are presented in Fig. [S] 
For all journals analysed the conclusions are common: the both hypothesis about pos- 
sible data approximations by log-normal ^ and by power-law with exponential cutoff 
(U) and a = 1 are almost equally good. Moreover, these conclusions are also right for 
the "Condensed Matter Physics" journal in spite of much smaller database with only 
262 records (see Fig. [G]). The observed data ffuctuations can be explained by relatively 
small statistics but such situation is usual for the majority of scientific journals. 

We could resume that it is hard to discriminate between the power-law or exponen- 
tial nature of the tail of the P{tw) distributions. But as is easy to see from obtained 
results the functional form of these distributions is the same: one explicit maximum 
and the long tail with several large values of tw The tail of the distributions built for 
several journals could be well approximated by functions with an exponential cut-off 
and leading power-law behavior t~^. Thus, we could consider the obtained form of 
probability distributions P{tw) as the typical one that can be used for scientometrical 
analysis of a given journals. 



2 Editorial process modeling 

We came to conclusion about the typical form of waiting time distributions P{tw) 
for scientific journals with normally working Editorial Boards based on the results 
described above. The origin of such functional form of P{tw) distributions is unknown. 
We can suppose that the contribution of human dynamics could be the reason of 
observed affinity of P{tw) to power-laws. The peer-reviewing stage (including the work 




8 



of authors and the communication process) is the human activity probably most similar 
to "natural" in comparison with the rest stages of editorial process. The other phases 
consist of different periodical tasks (Editorial Board meetings, uploads to web-site, etc.) 
or work with manuscripts in order of receiving (i.e., language and technical editing). 



I queue | ^ | queue ] 

>{I>-m>-^Irx^>o£>-|>^^^ 

typical time period ^[ ]]/]J/]|/7] ~ ] T ^{ [llTllTlilT "i 

[printedjssue] [printedjssue] 

Fig. 7 The schematic representation of modelled editorial work in scientific journals. 



To verify the hypothesis about the key role of peer-reviewing in the waiting times 
distributions P{tw) shaping we built the simple simulation model of editorial work 
in scientific journal omitting the peer-reviewing [IT]. Below we show the (expected) 
crucial difference of the waiting time distributions in such a model from the real data set 
analysed in the previous section. In the frame of this model, we consider the input flow 
of the submitted papers. Decision about manuscript acceptance is taken during regular 
meetings of the Editorial Board. All papers submitted before the meeting are considered 
(i.e. accepted or rejected) during this meeting. Generally, the Editorial Board meetings 
are held periodically and at least once before publishing of each issue. Let us defined the 
typical journal period T as a time interval (in days) between two consecutive journal 
issues. For example, since 1988 the journal "Physica A" has 24 issues per year and in 
this case T — 15 days. For "Condensed Matter Physics" journal (Fig. ^ T — 90 days 
since 1997 year. So, in our model the period between two consecutive meetings equals 
to the typical journal period T. Now the values of waiting time for each published 
paper could be calculated as the number of days between the moment of its receiving 
and the closest Editorial Board meeting. 

In the absence of peer-review all the received manuscripts are supposed to be 
considered by the Editorial Board. In this case the presence of some constrains is 
important (for example, the limited size of the issue). For simplicity all the papers in 
our model have an equal number of pages, so the maximal size of a single issue could 
be limited by number of paperj3- The value of parameter called the traffic intensity p 
(|2]) allows us to distinguish three regimes of work [5]: 

— the input flow of manuscripts is too slow - all the papers are accepted for publication 
but the issue is incomplete (subcritical regime, p < 1); 

— the number of received manuscripts is equal to the size of issue - all the papers are 
accepted for publication and the issue is complete (critical regime, p — 1); 

— the number of received manuscripts is larger than maximum number of papers in 
one issue - a part of papers are accepted and the other are waiting in the queue 
for the next issue (supercritical regime, p > 1). 

Obviously, the first regime of editorial work is not effective due to the incomplete- 
ness of issues. The third regime is also not realistic due to the endless queue of papers 
and therefore increasingly large values of waiting time. The second regime could be 



^ We specify the issue size equals 10 papers in our models. 
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considered as a perfect one since the issues are complete and the paper publishing 
without delays. But this regime is not stable and it could not be reached in practice 
because it is impossible to control the number of input manuscripts from different au- 
thors. It is more probable to provide the critical regime with X ~ fi, when the queue of 
papers periodically appears but it does not grow infinitely. So, the existence of limited 
queue is necessary to provide the completeness of issues being some kind of reserve. 

We start modeling with the simplest case when the number of input manuscripts is 
determined and equals to the size of issue (determined input flow). Besides, we set also 
the minimal possible size of printed issued. All the issues which are smaller than 80% 
from standard issue size are considered as incomplete. As manuscript submission times 
are uniformly distributed and all of them are accepted for publication, the resulting 
P{tw) distribution is uniform too, see Fig. [S] In this figure and further each point 
means the probability for the published paper to have waiting time in the interval 
tw G [t^..tl^^); the value of tw is measured in number of typical journal periods T. 
The results were obtained from the time-series simulations each of length 900 T. 



Fig. 8 The P(tw) distribution for the modeled periodical where the number of submitted 
manuscripts equals to the number of published papers (determined input flow). 



In practice the number of received manuscripts could not be always equal to the 
number of published papers. Traditionally, the input flow of tasks from an independent 
sources is modeled as the Poisson flow: 



Pk{t) = 



kl 



exp(-At), fc > 0, t>0, 



(5) 



here, fc is the number of input tasks during time interval t, X is the input flow intensity 
or the number of tasks which were received during the time unit. Poisson flow is also 
characterized by the exponentially distributed time intervals between two consecutive 
tasks. So, the modeling of Poisson flow mostly means generation of series of exponen- 
tially distributed random variables (for example, see 0). This approach is convenient 
for the cases when the execution law is also specified by the Poisson law. Then it is pos- 
sible to study different work regimes of the system controlling the values of intensities 
of input and execution laws. 

In the case of our model (omitting the peer-reviewing) we don't need to specify the 
execution law since all the papers could be accepted at one moment at the Editorial 
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Board meeting. So, the execution rate n could be expressed by the number of published 
papers while the number of received manuscripts could naturally define the arrival rate. 
Consequently, we need to simulate the input flow in a way allowing to control the input 
intensity by the number of received manuscripts. 

We generate the random value distributed by Poisson law (O with some value of A 
(input intensity) . In this case different number of manuscripts (which may be larger or 
smaller than A) could be received during the typical journal period T. Further, having 
the number of received papers we distribute them randomly over the T. If the number 
of manuscripts exceeds issue size, than a part of them goes to the queue or could be 
rejected. There are two ways of manuscript choosing from the queue (and also from 
the input flow) in our model: "FIFO" (flrst-in-first-out) and "RANDOM". We can 
suppose that the second way could be the simplest method to reflect the situation with 
the numerous continuous priorities. 




Fig. 9 The P(tw) distribution for the modeled periodical where input rate is comparably 
small: A = 3, = 10. Two scenarios of manuscript choosing from queue are used: (a) FIFO 
and (b) RANDOM. 



At first we modeled the case with non-determined input flow of papers and without 
any limitations of queue length. Modeling the input flow using ^ we can control its 
intensity changing the value of A. The modeling results are shown in Figs. 1^ and [HJl 
Here one can see the change of P{tw) distribution form according to an increase of A. 
If A is small then queue doesn't exist and waiting times are distributed more or less 
uniformly or forming the following "steps" . In this case the majority of published issues 
is incomplete (i.e., above 98% incomplete issues in the case of A = 3 and /i — 10). 

The queue appears when A is of the order of magnitude of fi. In this (critical) regime 
the queue length fluctuates around some value: the "excess" manuscripts from previous 
periods can flU the deflciency of papers for next periods. The P{tw) distribution in this 
case has two distinct "steps" (Fig. llOp presenting two categories of papers: which are 
printed in the flrst following issue and which are reserved for next periods in the queue. 
In our model such critical regime was reached at A = 7 and /i — 10. Then the number 
of incomplete issues descends down to 56%. 

To minimize part of incomplete issues it is enough to increase A still more. Our 
modeling results at A = 10 and fi = 10 are presented in Fig. 1111 This case corresponds 
to the unsteady state of the system. Despite the fact that number of incomplete issues 
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Fig. 10 The P(tw) distribution for the modeled periodical where input rate is approximately 
equal to the value of execution rate: X = 7, fi = 10. Two scenarios of manuscript choosing 
from queue are used: (a) FIFO and (b) RANDOM. 



is close to zero, such work regime is not efficient due to increasing queue length (one 
can see the queue length growing on the corresponding insets in Fig. [TT|). 




C - - - ^ " ' t 



Fig. 11 The P(tw) distribution for the modeled periodical with A = 10 and fi = 10. Two 
scenarios of manuscript choosing from queue are used: (a) FIFO (two cases) and (b) RANDOM. 
The queue length growing is shown on the corresponding insets. 

We can conclude that it is impossible to reach the optimal work regime for Edito- 
rial Board in scientific journal without any artificial constrains. So, the next steps in 
description of this process may involve additional constrains into the model. We limit 
the queue length to the issue size specified before. In our case normal issue volume 
equals 10, so the queue length is limited to 10 also. Thus, such work regime could 
be called "the issue in reserve". According to our rules the excess papers should be 
rejected using one of the scenarios: "LIFO" (last-in-first-out) and "RANDOM". 

As we can see from Fig. 1121 the obtained functional form of P{tw) distribution 
considerably differs from the experimental one. When the "FIFO" scenario is applied to 
select the manuscript from a queue the waiting time distribution is symmetric with one 
distinct maximum (Fig. 1121 (a) and (b)). In the case of "RANDOM" selection rule the 
waiting time distribution has decreasing form but also several "steps" corresponding to 
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each time period could be observed (Fig. [111(c) and (d)). Besides, the rate of decreasing 
observed is larger (with value of exponent close to 3) and more close to the exponential 
function. 

To summarize this section several aspects should be mentioned. First of all, it is 
important to note, that it is impossible to any real scientific journal to publish all 
the manuscripts received. So, every Editorial Board should choose the way or rule 
to eliminate the part of them: examining the execution of some obligatory technical 
requirements, applying peer-review mechanism, or following by other own internal cri- 
teria. The second important remark is that one can suppose that for different scenario 
it is possible to notice different picture of waiting time distribution for the published 
manuscripts. In this section we made the attempt to simulate most simple case of 
manuscripts arrangement in scientific journal, i.e. without any external peer review 
but using just simple rules. The obtained functional form of papers waiting time distri- 
butions for this scenario differs from the analogous results obtained for the real journals 
with peer-review mechanism. Of course, it is also interesting to simulate exactly the 
work with peer-reviewing but the way of natural human activity simulation is the 
subject of the discussions yet [3]. 
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Conclusions 

We have found that both log-normal and power-law function with exponential cutoff 
and exponent a = 1 can be the probable functions of waiting time distributions P(tw) 
for manuscripts in scientific journals. In fact, both log-normal and power-law functions 
predict exactly the same leading power behavior t^^, differing only in the functional 
form of the exponential correction TO^. In this sense, process considered here is gov- 
erned by similar probability distributions as other examples of human activities [HH 
131131015]. The observed data fluctuations can be explained by relatively small statis- 
tics but such situation is usual for the majority of real-world data bases. Thus, we 
consider the obtained form of probability distributions P{tw) as the typical one that 
can be used for scientometrical analysis of editions. The length of waiting times is 
an important characteristic of Editorial Board's activity. The publication delay effects 
journal rankings according to the impact factor [12] as well as personal citation rating 
of authors. 

The simple model of Editorial Board work was created to verify the hypothesis 
about the important role of peer-reviewing in the waiting times distributions P(tw) 
shaping. In general, the obtained results can be considered as the support of our pre- 
vious conclusions. 

In conclusion it is worth to mention some peculiarities of scientific editorial process 
that were not taken into consideration in our work. Here the main idea is connected 
with emphasizing of the meaning of natural human activity within the complex pro- 
cesses taking place in the Editorial Boards working. But these natural human activity 
patterns could be also affected by different situations or rules. For example, journal pol- 
icy could be very different concerning the strictness of time deadlines for peer-reviewing 
or concerning the manuscripts with one "positive" and one "negative" referee reports. 
In addition to peer review we can mention other processes as the examples of natu- 
ral human activity: technical work with raw manuscripts, communication between the 
editors and other participants of the editorial processes, etc. Moreover, one can distin- 
guish several human activity sub-processes within peer reviewing such as pure referee 
work with manuscript, communication processes between referee, editorial office and 
authors, and, eventually, the manuscript revision by authors according to the referee 
remarks. It is very complex problem to account all these aspects in one model. So, it 
could be the worth challenge for the future analysis. 
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